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The free fall of electric charges and dipoles, radial and freely falling into the 
Schwarzschild black hole event horizon, was considered. Inverse effect of electro- 
magnetic fields on the black hole is neglected. Dipole was considered as a point 
particle, so the deformation associated with exposure by tidal forces are neglected. 
According to the theorem, "the lack of hair" of black holes, multipole magnetic fields 
must be fully emitted by multipole fall into a black hole. The spectrum of electro- 
magnetic radiation power for these multipoles (monopole and dipole) was found. 
Differences were found in the spectra for different orientations of the falling dipole. 
A general method has been developed to find radiated electromagnetic multipole 
fields for the free falling multipoles into a black hole (including higher order mul- 
tipoles — quadrupoles, etc.). The electromagnetic spectrum can be compared with 
observational data from stellar mass and smaller black holes. 

1. INTRODUCTION 

It is known [1, 2j, black hole has no "Hair", so all multipole moments of the electromag- 
netic fields disappear as the system of charges close to the horizon of the Schwarzschild black 
hole. Field problem for static point charge has been solved by Linet In this problem it 
has been shown, in particular, that the field of a point charge close to the field of a charged 
black hole with the same charge, when the charge become close to the horizon. So all elec- 
tric and magnetic multipoles moments should be radiated, when the charge (or a system of 



* shatskiy@asc.rssi.ru 



2 



charges or currents) approach to the horizon of the Shwartzshild black hole. 

Loss in the energy is defined essentially by bremsstrahlung, when monopole's (unit charge) 
motion is accelerated. This is dipole radiation, since his power is inversely proportional to 
the c^, and radiated field components are inversely proportional to the (? . 

The existence of dipole radiation is unobvious in the case when the massive dipole fall into 
black hole, because both charges of dipole are moving and accelerated in the same direction, 
and the signs of the charges are opposite. However, increasing curvature of space leads to 
the radiation of the dipole type (see section |4]). 

Complete analogy can be made with quadrupole radiation of gravitational waves in the 
orbital motion of the two masses (whose power is inversely proportional to c^) in quadrupole 
approximation. In this case the signs of charges and masses are same, and velocities and 
accelerations are oppositely directed. 



2. THE LAW OF MOTION OF FREE FALLING PARTICLE 

Let's consider the radial free falling charge or electric dipole with mass m on the 
Schwarzschild black hole horizon. Schwarzschild metric has the form ^: 

Here Vg = 2M is the radius of the horizon of the Shwartzshild black hole. As you know, 
rapidly moving charge radiates. There is no radiation in the local neighborhood of the 
charge in their own co- moving frame ^. 

The acceleration of a massive particle ^ is defined by well-known expressions in a rigid 
system of reference at rest with respect to the black hole (see[4j, §87): 

^ Here the system of units is selected, in which c = 1 and G = 1 — the speed of hght and the gravitational 
constant. 

^ Co-moving freely falling frame of reference is not a rigid system, so there is always a relative acceleration 
between its different points, which can be neglected only in the locally small neighborhood of the charge 
(in Einstein's falling elevator). There is no radiation in this free falling "elevator" from free falling charge 
(in the approximation of small deformations of the reference frame in the "elevator" during his fall.) 

^ Mass of the falling particle assumed to be negligibly small compared with the mass of the black hole. 
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Here is component of 4-velocity and T^j are Christoffel symbols. We have for a 
Schwarzschild black hole (see 0], §100, 102) ^: 
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Here re(t) is the radius of the particle, tq is the radius at which the falling of the particle 
has been started and v{t) is a three-dimensional velocity of the particle in the Schwarzschild 
coordinates. 

The law of motion t(re) of the test particle is noticed in many papers ([2J, §2.4), which is 
free falling in the Schwarzschild field by radial trajectory from an infinite radius: 
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We need a similar law for a fall from a finite radius tq. 
The original formula is ^: 



(5) 
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So the initial conditions correspond to t{re = tq) = 0. 

Let's imagine that there are clocks showing Schwarzschild time t along the entire tra- 
jectory of the falling particle. Time At + Ati passed as long as the light signal from the 
clock (which the particle flies by) reaches the observer. The time interval At is the distance 
between the particle (when the particle flies past the clock) and its initial location (at the 
radius ro). The time interval Ati = const corresponds to the distance between the initial 
radius tq and a distant observer. 

According to [4j, §101 for period At we have 

^ro-rg' 



At = ro — Te + r„ In 



Tp - rn 



(7) 



^ Here and dot means differentiation with respect to timei. 

^ The similar formula is in the work |2j, but it is for dependence of falling particle's own time r on radius. 
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It is assumed here and below, that the distant observer is almost on the same line, as the 
particle. Dependence on a small angle 6 << n can be negleted (in linear approximation) 
between radius vectors of the observer and of the particle in the formula ([T]). 

Measurement of fall velocity of the particle by a distant observer must be described by 
time T(re) = t + At + Ati: 

_cos 

2 V^o 



T{re) = yr^{ro - r^){ro/rg - 1) + ^ 



+rg ■ In 
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+ const. (8) 
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Falling particle reach the horison in infinite time due to the logarithms in ([6]), ([T]) and (|8]). 
Let's introduce required notations: 

a = y^l - Tg/ve, b = \Jrg/re - Tg/rQ, Ci = ^1 - r^/ro . (9) 



(10) 



We have according to expressions (@ and (0): 

dT _l 1 _ h + ci 

dve fe 1 — fg/re a?h 

/?W-^ = -r^. (11) 

di + ci 

It is important to note, that /3 = for Tg = and Tg = (as well as for the quantity Tg). 
Physical meaning of the /3 is observable (by a distant observer) velocity of the falling 
particle. 

3. SPECTRAL DENSITY OF RADIATION 

We use different method, than in the work [5H9], to calculate the spectral density of the 
electromagnetic radiation. 

The spectral density of the electromagnetic radiation from the radial current (which flows 
along 6'e = 0) is given by ([4J, §66): 

d£^ = ^^-^ ■ ■ c/o ■ rfw , do = 271 sin 9 d9 . (12) 

Here E^(w) and H^(w) are spectral densities of the physical components of electric and 
magnetic fields, which are orthogonal to the direction of wave propagation (at infinity); 
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w = kt — the time component of the 4D null vector of photon kj, and the components kt 
and kg are integrals of motion for each of the emitted photon. 
General relativistic generalization of this formula is: 



dS iKW\9ee9tt\){KW\9ee9rr\) ^2 ^^3) 

47r^ 



The expression (13) at infinity (far from the black hole) can be rewritten using that com- 
ponents of the electric and magnetic fields are equal by magnitude in the electromagnetic 
wave: 



dS^ = - ■r^-do-dw = - -r^-do- dw. (14) 

Thus radiated Fourier components of the field F^* and F^"" at infinity must be inversely pro- 
portional to the square of the distance in order to satisfy the condition of energy conservation 
of electromagnetic waves. 

The total energy radiated to infinity can be obtained either by integrating over frequency 



of expression (14), or by integration of the Poynting vector over time. 

Spectral density tensor components F™" of the electromagnetic field and the contravariant 
components of the electromagnetic field F™" are associated with each other: 



exp [+iw)f:] dt, (15) 

exp [-int] dn. (16) 

Thus the spectral densities are obtained from the field components by the usual Fourier 
transform. 

Because we assume, that beginning of the fall of the particle is the time t = for a 
distant observer, in this moment the radiation field must be absent (-F^d ~ t < 0). 
This corresponds to half of the sum of the even and odd components of the radiation field, 
provided that the components we believe matching in absolute value. 

Crpmn I I rpmn I I rpmn I \ 
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Then the Fourier transformation of (15)-(16) for the radiation fields can be rewritten as: 



oo oo 

F"'" ^ (sin [wt] + cos [wt]) dt = j F""" cos 
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(17) 





In these formulas the sine corresponds to the odd component and the cosine corresponds to 
the even component of the radiation. 

4. RADIATION OF A CHARGE 

We use the second pair of general relativistic Maxwell equations to find components of 
the spectral density F^" of the electromagnetic field ^ (see [4J, §90): 



^^ (v^F™) = An^r , J"(r, t, ro) = V 5 [r - " m " ■ (19) 

— \/—QUl 

Here ^—g = sin 6* and b\r — Te(t)] is delta function with a singularity at the point, where 
it is the charge at time t. 



If r > To, the right side of equation (19) is always zero (4- vector current = 0), so we 



obtain following result by integrating of component j* of equation (19) over Q in [0, Q\. 



a 

j dr iV^F'') d9 + ^/^F^' = 



(20) 



We rewrite (20) in the linear approximation in the small angle ^ << tt as ^: 



(21) 



^ The sum over e is over all charges (for a dipole — over the first and the second charges) . 
^ Dependence of the field on the angle Q corresponds to the known relativistic expression of the work |10| 
or a, §67. 
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Here Ai is a 4-vector of the electromagnetic field of a falling charge in the stationary reference 
frame of a distant observer. Falling particle moves with a 4-speed [/* relatively distant 



observer. Therefore expression (11) helps us to find the form of U' 



Ti 8. 



W { —^=, , ^ , 0, } (22) 

Here and below the arrow indicates the limit as r/vg oo. 

We use two invariants to find the components At and Ar in the frame of distant observer : 
invi = (AiW) = At and inv2 = (AiAjg^^) = (At)"^. Here At is the only non-zero component 
of 4-vector Ai of the electromagnetic field in the free falling frame reference associated with 
the particle. The expression for At is known from work [3j: 



At = ^ 
rvp 



(r - M)(re - M) - M'^cosO ^ ^ 



v/(r - M)2 + (re - M)2 - M2 - 2(r - M) (r^ - M) cos 9 + M"^ cos^ 9 



(23) 



Here M = rg/2 is the mass of the Schwarzschild black hole. In the linear approximation 
respect to the angle 9 we obtain ^: 



T qejr - rg) 

At{e=Q) = r — 

r[r — re) r 



(24) 



This shows that the required asymptotic of capacity has form of the Coulomb field in the 
free falling co-moving frame of reference. This essential and logical result is a consequence 
of the local inertial frame of reference in the falling "Einstein's elevator." 

Then we obtain the required components with omitting all intermediate calculations: 

1 



At = UtAt ^ 
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Expression (21) can be written in the limit r/r^ — )■ oo as 



(25) 
(26) 



° We neglect the curvature of space-time at the point of a distant observer in the expression (22 1. 
^ Expression (23 1 depends on the small angle 6 only in the quadratic approximation. 

^'^ A member drAt has an asymptotic qe/r'^, so it decreases faster at infinity than a member dtA^ (which has 



asymptotics Qe/r). 
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Figure 1: Radial fall of a point charge q^, from radius vq = 5rg into the black hole. 

Curve with variable sign and greater amplitude is the spectral density (Fourier transform) of the 
radiation field of the argument (wr^) in units ^qed/r'^y 

Curve with fixed sign and lower amplitude is the square of F^, which is the power spectrum of 



electromagnetic radiation (14). 
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Here F^* is the tangential electric field of the radiation of a point charge when it is in the 
distance Vg from the black hole. It should be noted that this field is dipolar, since it is 
proportional to and inversely proportional to the square of the speed of light. 



We will get from (27) using the derivative with respect to time: 



/3 



2^3/2 



df3 _ _rga'^{h + ci)[a2ci - 2h^{h + ci)] 
d7g ~ 2cir2[(6 + ci)2 -a462]3/2 



We obtain by substituting (27) and (28) in the Fourier transform (17): 
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We obtain by replacing in (29) integration variable: dt = dre/re, according to Km for the 



spectral density of the radiation field: 
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eqe /■ w sin[wt - f ] /3 _ % f wci sin [wt(re) - fj 
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(30) 



We used here the integration by parts and we noted that /3 = in points = and = rg. 
The function t(re) is given by expression ([g]). 

Results obtained using (30) are shown in Figure [T| The spectral density of the radiation 
field of a charge is shown below. Here we use posibility of integration (29) in parts, according 



to the expression (26) and (30). 



■A. 



(31) 



Here Ai„ is the Fourier transform of the retarded vector potential for the field of the 
charge. It corresponds to the usual definition of the Fourier components of the radiation 
field (|E^| = wsin0|Aw| — see [4j, §66), and it is possible by ability of integration by parts 



in (29). But dipole field representation in the form (31) is not possible, as it will be shown 



in the next section. 



5. DIPOLE RADIATION 

The value of the dipole is defined as: do = qJq = const. We assume the dipole as point 
particle (/o <<< Vg) and neglect the deformation related to effects of tidal forces. We denote 
the radius of the dipole location: = (ri + r2)/2, then respectively replace — ?■ in all 
expressions for dipoles. 

5.1. Transverse orientation of the dipole 

First we choose the orientation of the dipole across the radius (the simplest case). Then 



Iq = r (1^9(1, so the expression (27) for the tangential electric field of the dipole radiation in 
this case will be rewritten as ^^: 



F^* =Af) — = ^ = ^ rga^{b + c,)[a'c,-2b\b + c,)] 



Formula ( 32 ) differs from the corresponding formula ( 27 ) only by coefficient and factor in the denomi- 



nator of the integrand. 
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Figure 2: Radial fall of a dipole from radius ro = 5rg on the black hole. 

For transversely oriented dipole (top): variable sign curve with greater amplitude is the spectral 
density (Fourier transform) of the radiation field F^^^ of the argument (wr^), in units of [(io/(rgr^)] . 
Fixed sign curve with lower amplitude is the square of , which is proportional to the power 



spectrum (14) of electromagnetic radiation. 
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Figure 3: The same thing for longitudinally oriented dipole is F^*^ , in units of 



The function F^^ij^^^) is the tangential electric field of the dipole radiation at the time of its 
location at a distance from the black hole (for transverse orientation). 
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So the spectral density of the radiation field of a dipole takes the form in the transverse 
orientation: 



cos [wt -j]^^_do f Tg cos[wt(rd) - f ] (6 + ci)[a'ci - 2h\h + ci)] 



4y26r3[(6 + ci)2-a462]3/2 



drd (33) 



Here the variable of integration was replaced with dt = drd/vd, as in the expression (30). 

It is important to note that the dipole radiation is independent of 6 in the linear approx- 
imation by the angle 6 << it in the case of the transverse orientation (since the radiation 
pattern of the dipole is located just in the transverse direction). 



5.2. Longitudinal orientation of of the dipole 

Now let's choose the orientation of the dipole along the radius. Length of the dipole can 
be calculated using following formula taking into account the curvature of space: 



Hence it follows 



/o^ h(t)-ri(t)]/y/l-r,/r. 



Ard = r2{t) - ri{t) = a ■ lo , Avq = ci ■ Iq. 



Here A denotes the change at constant time t, so differential must be zero: 



(34) 
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dt , dt ^ 
dt = — ■ Aro + — ■ Ard = 0. 



Expression (27) for a dipole can be written as: 
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Similarly to the case of the transverse orientation of the dipole the spectral density of the 
radiation field for the dipole in the longitudinal orientation takes the form: 
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The derivatives necessary for calculation of the quadrature (with regard to the expression 



(28)) are given in Appendix (see (39), (40) and (41)). 

Plots of the power spectrum for the transverse and longitudinal dipole orientations are 
shown in Figure [2] and Figure |3| respectively. 



6. DISCUSSION AND CONCLUSIONS 



As shown in previous sections, the radiation both falling monopole and falling dipole ((32) 



and (37)) are dipolar, because radiated field components are proportional to (28), that is 
inversely proportional to c^. 

The total energy of the radiation received by integration over time of the Poynting vector is 
£ = ^ J{F^^fr^dodt. It is equal ^ OmiGql/vg for the monopole and S'^ ^ O.OOOTrfg/r^ 
for the dipole, taking into account the angular distribution of power (see jlO] or §67). 
In addition, the energy of the radiation exceeds their rest energy for the falling electrons in 
the field of black holes with a mass less than It speaks about the infidelity of above 

calculations for radiation in the classical theory for primordial black holes (with a mass less 
than 10^^ g). In these cases it is necessary to use other theory (quantum gravity). 

Found spectrum of dipole radiation depends on the orientation of the dipole. 

Wavelength characteristic Am depends on the specific value of ro and on the orientation 
of the dipole at the maximum of the radiation. Am is approximately 20rg, when ro = Sr^. 
Therefore it's possible to observe this radiation near the maximum for a relatively small black 
hole (with a mass M <~ Mq). If you'll try to observe this radiation at shorter wavelengths, 
spectral power decreases approximately as w~'* in the local maximum . 

Since this radiation has a spectrum characteristic, it can be registered for the rare cases 
of falling magnetized planets (or asteroids) on black holes of stellar mass or even in rare 
cases, when pulsars fall into the black hole. 
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1. APPENDIX: necessary for the integration of quadrature (38) expressions 



We denote x = Tgjr^^ y = rgjr^ and give the necessary expressions: 
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4^2(1 _ ^)3/2^/^^ [l - x3 + x\y + 2) - 2xy + 2^ {\ - y){x - y) - y 
X Ix^^x - ?/ + (2a/1 - y - 7a/x - 

+ 2 (x^ - 2a; - 15) (V^^ + v^l - 2/) + (^9v/x^ - My^l - 
+ y\-x'' [y/^^ + 2v/r^) + x^ (2v/r^ - V^^) + X (29,/¥^ + 42 VT 
+ 37y^^ + 22v/r^|-xf29V^^+18v/r^) -7v/^^-2/r^ 
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X 2x^^x^ - (^17v/x^ + ISv^l -y^+x^ (Sy/x^ - 26^/1 - y^ 

- 4x^ (25y/¥^ + U^yi - + 16 {x^ + 4x - 3) y^ (^V^^ + v^T^) 

+ x^ (62^¥^ + 78y/l - + 2y^ (V^^ + v^l - y) - (25v/^^ + 29^1 - y) 

- x^ (^79v/^^ + 91v/l -y) + a; (337^^ + eiv^l -y) + 22v/^^ + 10^1 - y^ 
+ y |-x^ (y^^ + 2v/l-2/) + 12x' + 



+ 2xM50v/x^ + 71Vl -2/ -4xM3Vl -2/-26v/x^ - x ISl^x^ + 104 ^1 - 2/ 



+ 8y/x -y + \2^\ - 2/| + X [yJ^x -y- \^^\ - y) - 4^x - y 



(41) 



It is interesting that all these expressions are singular at the point r^^ = Tq, but their combi- 



nation in the formula (38), is finite. 



